Quantum information analysis of the phase diagram of the half-filled extended 

Hubbard model 
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We examine the phase diagram of the half-filled one-dimensional extended Hubbard model using 
quantum information entropies within the density-matrix renormalization group. It is well known 
that there is a charge-density-wave phase at large nearest-neighbor and small on-site Coloumb 
repulsion and a spin-density-wave at small nearest-neighbor and large on-site Coloumb repulsion. 
At intermediate Coulomb interaction strength, we find an additional narrow region of a bond-order 
phase between these two phases. The phase transition line for the transition out of the charge- 
density- wave phase changes from first-order at strong coupling to second-order in a parameter regime 
where all three phases are present. We present evidence that the additional phase-transition line 
between the spin-density-wave and bond-order phases is infinite order. While these results are 
in agreement with recent numerical work, our study provides an independent, unbiased means of 
determining the phase boundaries by using quantum information analysis, yields values for the 
location of some of the phase boundaries that differ from those previously found, and provides 
insight into the limitations of numerical methods in determining phase boundaries, especially those 
of infinite-order transitions. 

PACS numbers; 71.10.-w, 71.10.Fd, 71.10.Hf, 71.30. -|-h 



I. INTRODUCTION 

The one-dimensional Hubbard model is one of the 
earliest- and most-studied models of strongly correlated 
itinerant electronsii Its phase diagram at half filling was 
one of the first treated using quantum Monte Carlo for 
quantum lattice systems?^ one of the first modern numer- 
ical methods for these systems. Thus it is perhaps sur- 
prising that the details of the phase diagram of this model 
at half-filling have only recently been establishedi^i^i^ 

In this work, we determine the ground-state phase 
boundaries for the half-filled system using meth- 
ods based on quantum inf ormat ion i o d i d ^ calcu- 
lated using the density-matrix renormalization group 
(DMRG)J^iiiii^ These methods have a number of ad- 
vantages over other methods involving gaps, correlation 
functions, or order parameters . ^^d^ First, they involve 
only the properties of the ground state (or, in practice, 
a numerical approximation to the true ground state), 
and so can be calculated easily and accurately. Sec- 
ond, the behavior of the various quantum information 
entropies can be directly related to quantum critical be- 
havior, so that they provide unbiased indicators of quan- 
tum phase transitionsj^i^i^ Finally, they are quantities 
both intimately related to the fundamental approxima- 
tion of the DMRG and easily and naturally calculated 
within the DMRG algorithm ] ^^d^ Comparison with re- 
sults from other indicators of quantum phase transitions 
lends insight into the accuracy and limitations of both 
the indicators and the underlying numerical method. 

The one dimensional extended Hubbard model^S 
(EHM) is a relatively simple model which describes hop- 
ping of spin-i fermions between neighboring sites on a 
lattice. The fermions interact with each other through 



on-site and nearest-neighbor interactions, as described 
by the Hamiltonian 
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+ rii-iriii + V'^ niUi+i . (1) 

i=l 1=1 

Here cj^(ci(j) creates (annihilates) an electron with spin a 
on site i, — cl^Cia and rii = The parameter 

t denotes the strength of the inter-site hopping, U the 
strength of the on-site, and V the strength of the nearest- 
neighbor interaction. Since the nearest-neighbor and on- 
site interaction originate from the (repulsive) Coulomb 
potential in most materials, we consider only U,V > Q in 
this work. In general, we will consider finite lattices of 
size TV with open boundary conditions. 

Although the model is relatively simple, the phase di- 
agram of the half-filled system is more complicated than 
one might expect. When y = 0, the EHM reduces to the 
one-dimensional Hubbard model which is exactly solv- 
able using the Bethe ansatz»22, From this exact solution, 
it is known that the system is in a spin-density-wave 
(SDW) phase for all U > Q. This picture is reinforced 
by the effective model for strong coupling {U /t ^ 1), the 
Heisenberg model, which has a SDW ground state for 
antiferromagnetic spin exchange. One expects the SDW 
phase to remain the ground state for finite V , as long 
as U /V 3> 1. In the opposite limit, V/U ^ 1, one ex- 
pects the system to form a commensurate charge-density- 
wave (CDW) with an alternating spontaneously broken 
symmetry, i.e., sites with average occupation {ni) > 1 
alternate with sites with (n^) < 1. Such a picture is 
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obtained both in strong-coupling perturbation theory^! 
and in a weak-coupling renormalization group treatment 
(51-ology)!^ In fact, the phase boundary between these 
two phases is predicted to be at C/ — 2V, to lowest order, 
by both treatments. 

While a phase diagram containing only these two 
phases would not be particularly interesting in and of it- 
self, a prediction that a bond-order- wave phase (BOW), 
a phase characterized by a spontaneously broken sym- 
metry in which strong and weak bonds alternate, occurs 
near the CDW-SDW phase boundary at (approximately) 
U = 2V for intermediate coupling^^ has sparked new in- 
terest in the phase diagram. In addition, g-ology predicts 
the CDW-SDW phase transition to be continuous, while 
strong-coupling theory predicts it to be first-order, i.e., 
discontinuous. As a result, a number of papers on this 
issue have been published in the last few years, some con- 
taining contradictory resultsi^'^'^i^ Apparent agreement 
in the numerically determined ground-state phase dia- 
gram is present in two recent treatments.— However, in 
both of these calculations, the same quantities, the spin 
and charge exponents, the bond order parameter, and the 
spin and charge gaps, are used to determine the phase 
boundaries, and in the latter one consistency checks or 
realistic error estimates are not present. In view of the 
fact that both underlying numerical methods. Quantum 
Monte Carlo (QMC) and the DMRG yield results that 
are numerically exact for this system, it is not surprising 
that the subsequent similar analyses of the data yields 
essentially the same results. As we shall see, alternative 
ways of determining the phase boundaries yield more de- 
viation than is indicated by this agreement. 

In recent years, various types of quantum en- 
tropies have proven themselves to be a useful and 
accurate tool for determining the nature and crit- 
ical parameter values of quantum phase transitions 
(QPTs)^^'iaaiii2a6i25^^ Our goal is to combine 

these newly developed methods with accurate numerical 
calculations to independently determine the phase dia- 
gram. 

In this article, we calculate the one-site, two-site, and 
block entropies of the EHM using the DMRG and use 
these quantities to determine the nature and location of 
the phase boundaries in its phase diagram. The remain- 
der of this article is organized as follows: The concepts 
and methods used in our calculation will be discussed 
in Sec. [TTl with Sec. Ill Al defining the one-site, two-site, 
and block entropies and Sec. Ill Bl describing the DMRG 
method used. Results for these quantities, calculated us- 
ing the DMRG, will be presented in Sec. |TTI1 Sec. HITCl 
explores the limitations of using quantities such as order 
parameters to determine phase boundaries. Finally, the 
implications of our results as well as a comparison with 
previous work is contained in Sec. IIVI 



II. CONCEPTS AND METHODS 
A. Quantum Entropy 

Wu et al^ have argued that, quite generally, QPTs 
are signalled by a discontinuity in some measure of en- 
tanglement in the quantum system. One such measure 
is the concurrence;^ which has been utilized by a num- 
ber of authors^i^ i^°i^^i^^i'^°'^^i'^^ in their studies of spin 
models. The local measure of entanglement, the one-site 
entropy has been proposed by Zanardi^ and Gu et alX, 
or the negativity, by Vidal et alM- to identify QPTs. All 
of these quantities exhibit anomalies, i.e., either disconti- 
nuities in the quantity itself or in its higher derivatives or 
extrema, at points in parameter space that correspond to 
QPTs in exactly solvable models. What kind of anomaly 
occurs depends on the quantity and on the type of QPT; 
sometimes the anomaly appears only in some but not all 
of the quantum information entropies.— However, it has 
become clear from recent workiiii^ that, when the ap- 
propriate variant of the quantum information entropy is 
chosen, it can be used as a particularly convenient and 
accurate probe of the QPT. One example of this is the 
two-site entropy, a good probe of transitions to dimerized 
phases such as the BOW phase we study heieJ^ 

All of the quantum entropies discussed here have as 
their starting point the reduced density matrix of a bi- 
partite system. Assuming that a quantum mechanical 
system in a pure state 1 4') is divided into two parts, part 
A with complete orthonormal basis \i) and part B with 
complete orthonormal basis |j), the wave function can be 
expressed as 

= (2) 

The density matrix of the entire system is defined as 

P=l*>(*l, (3) 
and the reduced density matrix for subsystem A is 

PA = trB(|*)(^'|) (4a) 

= (4b) 

The von Neumann entropy for subsystem A is then 

Sa = -Tr pA^npA = -^Pa'^'a{pa) , (5) 

a 

where pa are the eigenvalues of the density matrix pA- 

It is also useful to define the density matrix in terms 
of an (unnormalized) projector which takes subsystem A 
from state to state \i), 
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For simple subsystems and an appropriate choice of ba- 
sis \i), the operator P^/ ^ can be expressed in terms 
of relatively simple operators, so that its matrix ele- 
ments can be directly calculated using the corresponding 
observablesj^^ this will be used in the following. 



1. One-site entropy 

In order to form the one-site entropy, which we will 
denote Se{l), the subsystem A is simply taken to be a 
particular single site £. Since Hamiltonian ([T]) contains 
no spin- flip processes, the reduced density matrix can be 
obtained directly in diagonal form if we take the spin oc- 
cupation basis (|0), I t), I i), I Ti)) as the basis states |i) 
in Eq. In this basis, we can directly obtain the eigen- 
values from the expectation values of the four operators^ 

(Ti \pe\ Ti) = {ne^riii) 
a \pe\ i) = ((1-nnWi) 
(T \pe\ T> = (n,T(l-nn)> 

The one-site entropy is relatively easy to calculate be- 
cause it requires only four local measurements on site 
£. While the one-site entropy is useful in some cases 
for characterizing first-order QPTs, it is typically not 
well suited for determining higher order QPTs because 
anomalies are often only discernible for large system sizes 
and are sometimes not present at aWj^ In particular, 
changes in intersite bond strength have no influence on 
it. Since the one-site entropy necessarily only depends 
on quantities that are localized on that one site (in fact, 
only on the average occupancy and on the average double 
occupancy), it cannot contain spatial information that is 
nonlocal. 



2. Two-site entropy 

It is therefore often useful to examine the von Neu- 
mann entropy associated with a larger subsystem. One 
convenient choice of subsystem A is that of two sites 
p and q, i.e., the two-site entropy Sp^q. In particular, 
we are interested in characterizing a BOW phase, which, 
for open boundary conditions, will have a broken bond- 
centered spatial symmetry. Thus, we are principally in- 
terested in the behavior of the two-site entropy for dif- 
ferent bonds, i.e., we take p and q to be pairs of nearest- 
neighbor sites. Therefore the two-site entropy will be 
denoted as 5^(2) in this paper, meaning Sp^p+i. 

The two-site entropy can be obtained by calculating 
and diagonalizing the reduced density matrix for the two 
sites. As for the one-site density matrix, its matrix ele- 
ments in the occupation number basis can be expressed 
straightforwardly in terms of expectation values localized 
to the two sites by considering the projector of Eq. ©• 



However, the result is necessarily somewhat more com- 
plicated than for the one-site density matrix. Since spin 
and particle number are conserved quantum numbers in 
Hamiltonian ([1]), the resulting reduced density matrix is 
block diagonal (rather than diagonal as for the one-site 
case) and has 26 independent matrix elements; for de- 
tails, see Ref. IsBI . Hence, 26 independent measurements, 
followed by the appropriate matrix diagonalization and 
the summation of Eq. ([5]), must be performed for every 
two-site entropy calculated. 

Note that calculating all N{N — 1) two-site entropies 
would be prohibitively expensive for large system sizes. 
While considering only nearest-neighbor bonds reduces 
this to — 1, it is usually sufficient to calculate the en- 
tropy of the two pairs of innermost sites, 5'Ar/2_i(2) and 
S'jv/2(2), i.e., two bonds chosen to be in the middle of the 
system to minimize boundary effects, to characterize a 
BOW phase. 



3. Block entropy 

The block entropy is also based on splitting the sys- 
tem into two subsystems A and B. However, subsystem 
A is now taken to contain the £ contiguous sites 1 to £ 
and subsystem B to contain the remainder, sites -I- 1 to 
N. (Note that other choices of sets of contiguous sites 
are also possible.) We then calculate the von Neumann 
entropy S{£) using Eq. © for £ G 1,...,N. Note that 
calculating the projector Pi'^i of Eq. ([6]) would be pro- 
hibitively complicated and expensive. Here, we instead 
calculate the density matrix directly from the wave func- 
tion using Eq. (|4cp . Fortunately, the wave function for 
all divisions of the system £ is readily available in the 
appropriate finite-system step of the DMRG algorithm, 
and, in fact, the corresponding von Neumann entropy 
is intimately related to the approximation made in the 

Unlike the two-site entropy, which has an finite upper 
bound (In 4 for Hubbard- like models), the block entropy 
generally grows as O (In N) for critical one-dimensional 
systems'^ - (but scales to a finite value for non-critical sys- 
tems). Although such a potential divergence would at 
first glance seem favorable for studying QPTs, the situa- 
tion is actually more complicated: boundary effects from 
open boundary conditions should have a stronger influ- 
ence on the block entropy than on a more locally deflned 
quantity such as the one-site or two-site entropy. It is 
therefore not clear which of these entropies can more ac- 
curately detect QPTs, but it seems sensible to expect 
that fast growing peaks are better detected by the block 
entropy, while non-diverging anomalies, such as disconti- 
nuities in derivatives, can be more precisely determined 
by the two-site entropy. 
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B. DMRG 

The DMRG calculations were carried out using the 
finite-system algorithnti^iiiii^ on systems with (mostly) 
open boundary conditions with from N = 32 to N = 512 
lattice sites. Open rather than periodic boundary con- 
ditions were used for two reasons: First, the DMRG 
is substantially more accurate for fixed system size and 
computational effort. Second, since open boundary con- 
ditions explicitly break the translational invariance, a 
corresponding spontaneously broken symmetry of the 
ground state, which, strictly speaking, can only occur 
in the thermodynamic limit, appears in the entropy pro- 
files of finite-sized systems. ^^''^^ This is, for example, the 
case for the BOW phase. 

Since only ground-state properties are required to cal- 
culate the von Neumann entropies, we need only calcu- 
late the ground-state wave function and the appropriate 
observables for the half-filled system. We have used dy- 
namic block-state selection (DBSS), which chooses the 
size of the Hilbert space retained in each truncation by 
keeping the block entropy of the discarded density-matrix 
eigenstates constant It is important to do this because 
the accuracy of the von Neumann entropy as calculated 
using the approximate DMRG wave function in Eq. (j4cp 
is directly related to the entropy threshold used. In our 
calculations we have used a threshold for the quantum 
information loss^ of x < 10~^°, which yields extremely 
precise results and requires a maximum of approximately 
m — 3000 block states to be retained. In general, we esti- 
mate that the errors due to the truncation in the DMRG 
calculations are negligible in comparison to the uncer- 
tainties arising from the finite-size scaling. 

C. Limitations of other quantities 

Naively, the most straightforward way to determine the 
existence and extent of a BOW phase would probably be 
to examine the bond order parameter, defined as 

1 

■B = ^ XI (^1+1'^'=''^ + 4aC^ + la) ■ (7) 
i— l.cr 

Whether the bond order parameter is finite or vanishing 
in the thermodynamic limit would determine whether a 
given point in parameter space is ordered or not, and a 
grid of such points can be used to determine the phase 
boundaries. 

Unfortunately, there are two major problems with this 
strategy. First, the transition between the BOW and 
SDW phases is expected to be infinite-order, and we in- 
deed find behavior characteristic of an infinite-order tran- 
sition. This means that the extrapolated bond order pa- 
rameter tends to zero exponentially as the BOW-to-SDW 
transition is approached. Second, while it is know that 
the bond order parameter is linear in 1/N in the CDW 
phase, where N is the system size, and proportional to 



1/VN in the SDW phase, the analytic form of the finite- 
size scaling in the BOW phase is not known and changes 
nature as the transition is approached. We will examine 
this issue in more detail in Sec. IIII CI 

Unfortunately, similar difficulties arise in other quan- 
tities that are typically used to determine critical pa- 
rameters. In particular, the charge gap goes to zero in 
a clear manner at the CDW-BOW phase transition, but 
does not exhibit any anomalies at the BOW-SDW transi- 
tion. The spin gap, on the other hand, does go to zero at 
the BOW-SDW transition, but, like the bond order pa- 
rameter, with an exponential dependence. Therefore, it 
is also not well-suited for exactly determining the phase 
boundary. 

III. RESULTS 

We consider first the behavior of the two-site and block 
entropies. In order to map out the phase boundaries, we 
have swept V from V = 2toV = 5.5 in steps of 0.5 for a 
number of values of U. For all sweeps, we find that two 
peaks develop in both the two-site entropy, Fig.[Tl and in 
the block entropy. Fig. [21 for sufficiently large systems. 
Note that the second peak in the block entropy can only 
be seen at system sizes of = 96 and larger. This slow 
size dependence might be the reason why this phase was 
not seen by Deng et. alr^. 
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FIG. 1: (Color online) Two-site entropy S'jv/2+i(2) for V/t = 
3, plotted as a function of U/t. System sizes range from = 
64 to iV = 512 sites. 

We interpret the peak at lower U as marking the 
CDW-BOW phase transition at the corresponding sys- 
tem size, and the peak at larger U as indicating the 
BOW-SDW phase transition. The differing shapes of the 
two peaks are consistent with the picture that the CDW- 
BOW transition is first order in this parameter regime;^ 
while the BOW-SDW transition shows characteristics of 
an infinite-order transition. 

The interpretation that the intervening phase is a 
BOW phase is supported by the behavior of the two-site 
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FIG. 2: (Color online) Block entropy at the center of the 
chain S{N/2 + 1) for V/t = 3 and system sizes from A'^ = 64 
to TV = 512 sites. 



entropy for two adjacent (odd and even) bonds in the 
center of the lattice, plotted for the largest system size, 
TV = 512, in Fig. [TJ The dimerization entropy is given 
by Dg = 5'jv/2+i(2) — S'Ar/2(2), which is the difference be- 
tween the two TV = 512 curves. It is clear that this differ- 
ence reaches a marked maximum between the two peaks 
associated with the phase transitions. Finite-size extrap- 
olation (not shown) indicates that remains finite in 
the thermodynamic limit in the intermediate phase. 

The positions of the two peaks in the two-site or block 
entropy can then be extrapolated to the thermodynamic 
limit. While the functional form of this extrapolation 
is not exactly known, using a fourth-order polynomial 
yields stable results, with a rapid falloff in coefficient size 
for higher orders. Therefore, the behavior is predomi- 
nantly linear. This can be seen in Fig. [3) 



A. The CDW-BOW transition 



A comparison of Figs. [T] and [2] clearly shows that the 
peaks in the block entropy at lower U are sharper and 
higher than those in the two-site entropy. In addition, 
the position of the peaks in the block entropy match the 
jump in the one-site entropy (see Fig. |4]) better than the 
position of the peak in the two-site entropy. Therefore, 
we conclude that the block entropy is, in general, a bet- 
ter indicator of the position of the CDW-BOW transi- 
tion than the two-site entropy. As can be seen in Fig. [31 
the fits to the two-site and the block entropies (the two 
lower curves) match almost exactly in the thermody- 
namic limit, so that this issue is virtually irrelevant here. 

A more difficult issue is to determine where nature 
of the CDW-BOW transition changes from first-order 
to continuous. This can best be investigated using the 
finite-size extrapolation of the one-site entropy. In Fig. IH 
we show the one-site entropy extrapolated to the thermo- 
dynamic limit using a polynomial of cubic order. It is ap- 
parent that the entropy has a jump at the transition point 
when V/t > 4, a clear indication of a first-order transi- 
tion. At approximately V/t = 3, the transition is close to 
becoming continuous in that a jump is no longer present. 
For smaller values of V/t (not shown), it is clearly con- 
tinuous. Therefore, we conclude that the first-order-to- 
continuous bicritical point must occur somewhere near, 
but below V = 3. Note that it is difficult to determine 
the location of this point with more accuracy because 
one would have to determine whether or not an increas- 
ingly small jump is present in the finite-size extrapolated 
data as the bicritical point is approached on a sufficiently 
finite grid. 
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FIG. 3; (Color online) Finite-size extrapolation of the peaks 
in Figs. [T] and [2] to the thermodynamic limit, 1/L — » 0, using 
fourth order polynomials in 1/L. "2s" labels the two-site 
entropy and "B" the block entropy. 
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FIG. 4: (Color online) Extrapolated one-site entropy 5'jv/2(l) 
plotted as a function of U/t for various values of V/t. 
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B. The BOW-SDW transition 

The BOW-SDW transition is believed to be infinite- 
order— and is therefore much harder to characterize. 
While we were able to obtain a fairly good estimate for 
the position of the CDW-BOW transition just by exam- 
ining the bond order parameter, this is not possible for 
the BOW-SDW transition, as discussed in Sec. IIII CI As 
can be seen in Figs. [1] and [21 the maxima for the BOW- 
SDW transition (the peaks at higher U /t) are much 
broader. Therefore, small errors in the numerical calcu- 
lations would have a bigger influence on the result than 
for the case of the CDW-BOW transition. Due to this, 
we have carried out very high precision calculations, as 
already described in Sec. IIIBI 

Results for the phase boundaries are shown in Fig. O 
plotted in the tilted V-U phase, i.e., with axes 2V/U and 
U , so that the transition region is discernible. Included 
are data from Ref . |3|, in which the phase boundaries were 
determined from the spin and charge exponents calcu- 
lated using QMC methods. There is generally very good 
agreement for the location of the CDW-BOW transition, 
with some deviation of the results from the two-site en- 
tropy at smaller U /t values. As we have argued above, 
the block entropy is a better indicator of the position of 
this transition because the peak is better developed and 
grows more rapidly with system size. 
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FIG. 5: (Color online) The phase diagram in the tilted U-V- 
plane, showing phase boundaries determined using the two- 
site entropy (2s), the block entropy (B), and including results 
from Ref. ^ determined from QMC calculations. 



in localizing the broad peak in the entropies correspond- 
ing to the BOW-SDW transition. In addition, due to 
the strong increase of the block entropy with system size 
below U/t « 4, we cannot treat systems of more than 
a few hundred lattice sites for fixed x = 10"^", leading 
to increased uncertainty in the finite-size extrapolation. 
For U/t < 3 the second peak in the entropy functions 
develop only for A'^ > 96, quite severely limiting the ex- 
trapolation to infinite system size. We therefore display 
the phase diagram only for U/t > 3 in Fig. [S] 



C. Bond-order-parameter results 



As already mentioned in Sec. Ill Cl we have carried out 
calculations for the bond order parameter to very high 
precision. The resulting data are sufficiently accurate 
so that they can be regarded as essentially exact for a 
particular size for fitting purposes. We have carried out 
the extrapolation to the infinite-system limit by fitting to 
three different functions: a polynomial in 1/N, a poly- 
nomial in 1 / -s/a", and a power law of the form 1 /N" . 
The extrapolated data are shown in Fig. [51 In the CDW 
phase, to the left of the transition indicated by vertical 
lines, the fit to a polynomial in 1/A^ gives the best result, 
yielding the expected value, zero, to the best accuracy. 
In the SDW phase, for U/t > 8, the polynomial fit in 
1/VN and the fit to a power law work better, yielding 
the expected value of zero to within reasonable accuracy. 
This is indicative of a scaling whose dominant term falls 
off more slowly than 1/A^. In the intermediate region, 
i.e., in the BOW phase, the results differ significantly, 
with both the fit to powers of 1/ ^/N and the power-law 
fit extrapolating to spurious negative values. In addition, 
the power-law fit is clearly unstable in the BOW region. 
While the fit to a polynomial in 1/A^ seems to be more 
stable, it clearly overestimates the bond order parameter 
significantly in both the SDW phase and in most of the 
BOW region. 

Therefore, we conclude that the behavior of the bond 
order parameter can be used to confirm the position of 
the CDW-BOW phase transition but is notoriously un- 
reliable for determining the location of the BOW-SDW 
phase transition. 



IV. PHASE DIAGRAM AND DISCUSSION 



For the BOW-SDW transition, the two-site and the 
block entropies coincide perfectly upon scaling for higher 
U and V values, where only a first-order transition is 
present, as can be seen in Fig. [5l However, there is a 
discrepancy in the position of the peaks in the infinite- 
system extrapolations of the two-site and the block en- 
tropies (see Fig. [3]) for smaller values of U and V. This 
occurs for the entire range of parameter values in which 
a BOW phase and, therefore, an infinite-order transition 
is present. This is probably partially due to uncertainty 



We now summarize the current state of knowledge of 
the phase diagram of the half-filled extended Hubbard 
model and discuss open issues and uncertainties. The 
overall phase diagram is relatively well understood; our 
results are depicted in Fig. [71 The presence of the SDW 
phase at F < 2C/ and the CDW phase at F > 2?/ have 
been long understood, as well as the fact that the transi- 
tion occurs a,t V K 2U. The picture of there being a sin- 
gle first-order transition line at strong U and is also 
well-established. Our work lends support to a picture 
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FIG. 6: (Color online) Bond order parameter B extrapolated 
to the thermodynamic limit using fits to three different func- 
tions: fit to a third-order polynomial in 1/A'' (Poly), fit to 
a third-order polynomial in l/\fN (Sqrt-Poly), and fit to a 
power law l/N" (Power-law) with a a fitting parameter, plot- 
ted as a function of U/t. The vertical lines indicate the tran- 
sition point determined using the two-site entropy (2s) and 
the block entropy (B). 



in which an intermediate BOW phase is present between 
the CDW and SDW phases for intermediate to small U 
and V\ our results indicate that this phase is present for 
V/t < 5. At this point, we find that the first-order CDW- 
SDW transition line bifurcates into a first-order CDW- 
BOW transition line and an infinite-order BOW-SDW 
transition line. The CDW-BOW transition line remains 
first-order at a bicritical point at somewhat smaller V/t, 
below which it becomes continuous, presumably second- 
order. These results are in reasonable agreement with 
the results of Refs. i and i. We therefore regard these 
features of the phase diagram as being well-established. 

We now discuss details of the phase transition more 
quantitatively. As we have seen in Fig. [5l there is not 
much uncertainty in the position of the CDW-BOW 
phase transition. The remaining interesting question for 
this transition is the location of the bicritical point, i.e., 
exactly where the phase-transition line goes from be- 
ing first-order to being continuous (presumably second- 
order). However, as we have pointed out in Sec. IIII Al 
entropy measurements can only roughly determine that 
this point occurs at around V/t = 3 and are not an ideal 
measurement to locate it more accurately. While other 
authors have obtained putatively more accurate values 
for the location of this bicritical point j^i^ we point out 
that the inaccuracies in our method reflect intrinsic limi- 
tations of the numerical methods, which stem both from 
the DMRG truncation error as well as from the limita- 
tions of working with finite systems. 

The exact position of the BOW-SDW phase-transition 
line is also somewhat uncertain. There is a small, but 
significant, discrepancy in our calculations between the 
values obtained from the two-site-entropy and those ob- 
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tained from the block-entropy; however, both of these 
extrapolated values seem to converge smoothly to the 
same line at the tricritical point. The most likely ex- 
planation for the discrepancy in the extrapolations lies 
in the finite-size extrapolation, i.e., more precise results 
could be obtained if larger system sizes could be treated. 
The deviations between the two values for the transition 
line can be taken as a rough estimate of the uncertainty 
in the position of the line. In addition, both of these 
values deviate from those of Refs. d andU (see Fig. E]), a 
deviation to larger values of V, i.e., to a narrower BOW 
phase, in both cases. 
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FIG. 7: (Color online) The phase diagram in the U-V-Plane 
obtained from our calculations. On this scale, uncertainties 
in the position of the transition lines are smaller or about the 
size of the symbols. The phase in the narrow region between 
the two transition lines is BOW. 

In summary, the DMRG method, coupled with the use 
of single-site, two-site, and block entropies, is a power- 
ful, relatively unbiased method to determine subtle prop- 
erties of phase diagrams such as that of the extended 
Hubbard model at half filling. However, the limitations 
of the numerical results reflect the intrinsic limitations 
of the method and the problem studied. In particular, 
three aspects of the phase diagram studied here remain 
difficult to pin down numerically: the exact position at 
which the CDW-SDW fine bifurcates into CDW-BOW 
and BOW-CDW transition lines, the position of the bi- 
critical point at which the CDW-BOW transition goes 
from first to second order, and the exact position of the 
infinite-order BOW-SDW transition. In addition, our 
calculations make clear that quantitative determination 
of the transition lines become very difficult in the region 
of small U and V . 
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